Diamagnetism and Cooper pairing above Tc in cuprates 



o 
o 

(N 
o 

Q 



o 
o 

Oh: 



o 



> 
m 

oo 

o 

O 



X 



Lu Li^, Yayu Wang^'^, Seiki Komiya'^, Shimpei Ono^, Yoichi Ando'^-'*, G. D. Gu^ and N. P. Ong^ 
^Department of Physics, Princeton University, Princeton, NJ 08544 
^Department of Physics, Tsinghua University, Beijing, China 
^Central Research Institute of Electric Power Industry, Komae, Tokyo 201-8511, Japan 
"^Institute of Scientific and Industrial Research, Osaka University, Iharaki, Osaka 567-0047, Japan 
^ Brookhaven National Laboratories, Upton, NY 11973 
(Dated: December 16, 2009) 

In the cuprate superconductors, Nernst and torque magnetization experiments have provided 
evidence that the disappearance of the Meissner effect at Tc is caused by the loss of long-range 
phase coherence, rather than the vanishing of the pair condensate. Here we report a series of torque 
magnetization measurements on single crystals of La2_xSrxCu04 (LSCO), Bi2Sr2_yLayCu06 (Bi 
2201), Bi2Sr2CaCu208+i (Bi 2212) and optimal YBa2Cu307. Some of the measurements were taken 
to fields as high as 45 T. Focusing on the magnetization above Tc, we show that the diamagnetic 
term Md appears at an onset temperature T^^^^t high above Tc. We construct the phase diagram 
of both LSCO and Bi 2201 and show that T^naet agrees with the onset temperature of the vortex 
Nernst signal To„sef Our results provide thermodynamic evidence against a recent proposal that the 
high-temperature Nernst signal in LSCO arises from a quasiparticle contribution in a charge-ordered 
state. 

PACS numbers: 74.25.Dw, 74.25. Ha, 74.72.Hs 



I. INTRODUCTION 



A series of experiments utilizing the Nernst Effect 
P, H, H, 0, Q has demonstrated that an enhanced Nernst 
signal is observed in hole-doped cuprates at temperatures 
T significantly above the superconducting transition tem- 
perature Tc- The high-T Nernst region was identified as 
a continuous extension of the vortex liquid state. In this 
strongly fluctuating vortex-liquid state, the large Nernst 
signal arises from phase slippage caused by singular phase 
fluctuations of the pair condensate @, [1] • In the phase- 
disordering scenario, the unbinding of vortex-antivortex 
pairs (in zero applied H) leads to the loss of long-range 
phase coherence at Tc [6|. The condensate is incapable 
of displaying long-range supercurrent response. Hence, 
even in weak H, there is no Meissner effect above Tc de- 
spite the survival of the pair condensate. Nevertheless, 
the persistent short-range phase stiffness supports vortic- 
ity and produces a large, strongly T-dependent Nernst 
signal in the presence of a temperature gradient — VT 
and an applied magnetic field H. The Nernst effect above 
Tc has also been investigated in Refs. 0, H, [1] 

Subsequently, thermodynamic evidence for the pair 
condensate above Tc was obtained by torque magnetom- 
etry, which is a very sensitive probe of diamagnetism in 
the cuprates [13, El- A large diamagnetic response, that 
is non-linear in H and grows strongly with decreasing 
T, is specific to the Cooper pair condensate. Conse- 
quently, the diamagnetism results present clear evidence 
for survival of the pair condensate - with sharply reduced 
phase stiffness - to temperatures high above Tc. This 
complements the transport evidence from the Nernst ex- 
periments. To date, the high-resolution torque measure- 
ments above Tc have been reported in the bilayer cuprate 
BisSrsCaCusOg+a M MM, MM and in lightly- 



doped La2-xSrxCu04 at low T [15]. For earlier mea- 
surements of diamagnetism above Tc, see Ref. [l^ (on 
Bi 2212) and Ref. [l3 (LSCO). 

We report further torque magnetization ex- 
periments on La2-xSrxCu04, Bi2Sr2_yLayCu06, 
Bi2Sr2CaCu208+5 and YBa2Cu307, which show that, 
when a positive Nernst signal appears below the onset 
temperature Tonset, it is accompanied by a large dia- 
magnetic signal that grows steeply with decreasing T. 
Extending the torque measurements to intense fields (33 
T to 45 T), we also show that the curves of M vs. H are 
nonlinear with a profile characteristic of vortex liquid 
response, even at elevated T. 

Recently, the fiuctuating vortex-liquid interpretation 
of the Nernst effect has been challenged by Cyr-Choiniere 
et al. [Til, who carried out Nernst measurements on a 
series of Nd- and Eu-doped LSCO cuprates, in which 
charge ordering associated with stripe formation is known 
to occur at a charge ordering temperature Tco- By 
a qualitative comparison of the Nernst coefficient v in 
Nd-LSCO and Eu-LSCO, Cyr-Choiniere et al. proposed 
that, in pure LSCO, the high-temperature Nernst sig- 
nal Cn arises from small quasiparticle pockets as a re- 
sult of Fermi-Surface (FS) rearrangement, rather than 
from phase-slippage in the pair condensate. We discuss 
the problem of separating the vortex Nernst term from 
quasiparticle contributions, and the key role that dia- 
magnetism plays in this task. In pure LSCO, we explain 
how the magnetization results reported here pose serious 
difficulties for this hypothesis. 

We adopt the abbreviations LSCO 09, LSCO 12 and 
LSCO 17 for La2-xSrxCu04 with the Sr content x = 0.09, 
0.12 and 0.17, respectively. We use Bi 2201 and Bi 2212 
to stand for Bi2Sr2_yLayCu06 and Bi2Sr2CaCu208+5, 
respectively, and YBCO for YBa2Cu307. The terms un- 
derdoped, optimally-doped and overdoped are abbrevi- 
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ated as UD, OPT and OV, respectively. 



II. TORQUE MAGNETOMETRY 

The torque magnetization was measured with the sam- 
ple glued to the tip of a thin cantilever, with H applied at 
a tilt angle — 10-15° to the crystal c-axis. The deflec- 
tion ip of the cantilever by the torque was detected capaci- 
tively. The cantilever typically can resolve changes in the 
magnetic moment of 5m ~ 10~^ emu. Because of the 2D 
electronic dispersion in cuprates, the diamagnetic orbital 
currents are largely confined to the a-h plane [13, E^. 
This makes torque magnetometry well-suited for detect- 
ing weak, incipient diamagnetism in cuprate crystals. We 
generally report the raw data as the effective torque mag- 
netization Meff defined as M^^ff = t/ fiQHV sinO, where 
T is the torque signal, /ip the vacuum permeability, and 
V the sample volume. {M^ff includes all contributions 
to the observed torque signal.) 
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FIG. 1; (color online) Magnetization curves of Sample LSCO 
09 with Sr content x — 0.09 and transition temperature Tc = 
24 K measured in magnetic fields H up to 14 T. (A) The (to- 
tal) effective magnetization Me// vs. H at temperatures 4.2 
< T < 200 K. Below Tc, the Me// - H curves are hysteretic 
when H lies below the melting field Hm{T), which is exper- 
imentally defined as the field at which the up-sweep branch 
deviates from the down-sweep branch {Hm is indicated by an 
arrow for the curve at 7 K). Above Hm, the curves become 
reversible. (B) The diamagnetic magnetization Md vs. H at 
temperatures 4.2 < T < 30 K. (C) Curves of Ma vs. H at 
22 K < T < 80 K displayed in expanded scale. In LSCO 09, 
the diamagnetic signal persists to more than 60 K above Tc- 
In (B) and (C), the bold curve is measured at the separatrix 
temperature Ts = 22 K. 

To describe our analysis, we first examine the curves of 
Meff VS. H shown in Fig[TK for Sample LSCO 09 with 
Tc = 24 K. At temperatures T > 100 K, Mejf is strictly 
linear in H and paramagnetic in sign. This reflects the 
dominance of the anisotropic Van Vleck paramagnetic 



susceptibility Axp, which has a weak T dependence given 
by Axp = A + BT, with A > BT > at 200 K. (The T- 
dependence of Me// will be shown later in Fig. [9l) The 
reason for measuring a dense set of curves in the high-T 
interval 100-200 K will emerge when we discuss extrac- 
tion of the onset temperature for diamagnetic response 
(Sec. IVI|) . Below 100 K, Me// begins to display a weak 
diamagnetic contribution that rapidly increases in mag- 
nitude as T decreases. The temperature at which the 
diamagnetic contribution (referred to as Md hereafter) 
appears is identified as the onset temperature T^nset- 

We note that M^ff becomes increasingly nonlinear in 
H as T decreases from 60 K to Tc (24 K). Below Tc, 
the diamagnetic term becomes so dominant in magni- 
tude that Me// is forced to large, negative values despite 
the positive contribution of the paramagnetic Van Vleck 
term. At T < Tc, the Afe// — H curves at low fields 
display strong hystereses between up-sweep and down- 
sweep branches because of strong pinning in the vortex- 
solid state. In the plotted curves, this is seen as a fork 
(we suppress the full hysteretic curves in both branches 
for clarity). The field at which the down-sweep branch 
deviates from the up-sweep branch is experimentally de- 
fined as the melting field Hm{T) of the vortex solid (indi- 
cated by arrow for the curve at 7 K). As T is raised, Hm 
decreases rapidly, reaching zero at Tc. A detailed inves- 
tigation of Hm{T) in hghtly doped LSCO is reported in 
Ref. [IHl • The spin contribution to Me/ / , which becomes 
important below 2 K in samples with x < 0.06, is also 
treated in detail in Ref. [ISi] . 

We assume that the paramagnetic background AxpH 
follows the trend that is seen at T > Tonset- Hence the 
diamagnetic term M^ is related to the observed torque 
magnetization Mcff by 

Mcff{H) = Mci+Axp{T)H - MciiH) + {A+BT)H. (1) 

Hereafter, we subtract the background Van Vleck term 
{A + BT)H from Me// and plot the diamagnetic term 
Md vs. H (except stated otherwise). 



III. LSCO 

Carrying out the background subtraction for LSCO 09, 
we obtain the Md vs. H curves. They are displayed at 
selected T in Fig. [Tp (4.2-30 K) and Fig. [Tp (20-80 K). 
As shown in Panel B, M^ is nonlinear in H over a broad 
temperature interval. The curve at 22 K (bold curve) 
displays a characteristic fiat profile in low fields (0.5 T 
< H < 5 T). We identify this temperature as the "sepa- 
ratrix" temperature Ts [1(|. Below Tg, Md takes on very 
large, negative values at small H . As H increases, Md 
displays an initially steep logarithmic increase, followed 
by a slower approach towards zero as H approaches the 
upper critical field Hc2- The low-field curvature of the 
Md vs. H curve changes from negative below Tg to pos- 
itive above T,. To emphasize the high temperature dia- 
magnetic response, panel C displays the Md — H oi Sam- 
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pie LSCO 09 at T > Ts in expanded scale. For T > T^, 
the curves remain diamagnetic, displaying pronounced 
nonlinearity vs. H. We regard the nonlinear diamagnetic 
response above Tc as clear evidence for the presence of 
local supercurrents as well as finite pair amplitude in the 
pseudogap state. 




10 20 30 

^0 H ( T ) 



FIG. 2: (color online) Magnetization curves Md vs. _ff in 
LSCO 09 measured in intense fields up to 33 T. Tc lies between 
the curves at 22 K (bold) and 25 K. 
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As in the superconducting state, the diamagnetic sig- 
nals above can be suppressed by an intense magnetic 
field. To get a sense of how large this field scale is, we 
extended the torque measurements on LSCO 09 to 33 T. 
After the background subtraction described above, the 
resulting Md — H curves are isolated and plotted in Fig. 
[21 For T > Tc, these nonlinear — H curves display a 
broad minimum. The characteristic field Hmm locating 
the minimum increases rapidly with T (from 8 T at 25 
K to 33 T at 40 K). A rough measure of the field-scale 
needed to observe the inherent nonlinear response above 
Tc is given by Hmin- 

We emphasize that, above 50 K, Md is seemingly linear 
in H. However, this is simply because Hmm has now 
moved outside the experimental window. It is incorrect 
to attribute - using a limited field range H <15 T - 
this putative linear behavior to a different mechanism 
(e.g. quasiparticles with Landau diamagnetism) . Above 
Tc, Md is inherently non-linear in M up to the onset 
temperature T^l^^t^ but one needs progressively higher 
fields to see the non-linearity as T increases. 

The diamagnetic Md — H curves of Sample LSCO 17 {x 
= 0.17, Tc = 38 K) and Sample LSCO 12 [x = 0.12, Tc = 
27 K) measured in H up to 45 T are plotted in Fig. 
and Fig. respectively. In these higher-doped crystals, 
the pattern of the Md — H curves is broadly similar to 
that in LSCO 09, and the foregoing discussion apphes 
to the diamagnetic curves. The higher field accessed (45 
T) in Fig. [3J'V, B confirms the intrinsic non-linearity of 
Md{H). In Panel A, Md is manifestly nonlinear at all T 
up to 70 K, even though it seems linear when the field 
range is restricted to 7J < 10 T. 

The magnitude of Md in the high-field curves reveals 
an interesting difference (for T < Tc) between LSCO 17 
and LSCO 12. The monotonic decrease of \Md\ with 
H provides an estimate of the upper critical field scale 
Hc2 (by extrapolating Md — > 0). Comparing LSCO 17 
with LSCO 12, we see that the scale of \Md\ above 20 
T is 2-3 times larger in the former at the same H and 
T. However, the stronger curvature of Md vs. H in 
LSCO 12 implies that the decay of Md is more gradual, 
so that its Hc2 is actually higher than that in LSCO 17. 
In the related cuprate La2-xBaxCu04, stripe formation 
at a; = i drives to 4 K 20]. The "dip" in the Tc dome 
> in the phase diagram of LSCO suggests that fluctuating 
stripes may also exist at a; = |. If this is true, the Md 
results suggest that fluctuating stripes reduce the overall 
pair condensate strength (compared with x = 0.17) but 
allows it to survive to slightly larger fields. 



IV. BISMUTH-BASED CUPRATES 



FIG. 3: (color online) Comparison of magnetization curves 
Md vs. H in LSCO 17 and LSCO 12 measured in intense 
fields up to 45 T. Panel A shows curves for Sample LSCO 17 
(x — 0.17 and Tc = 38 K). Panel B displays curves of Sample 
LSCO 12 {x = 0.12 and Tc = 27 K). 



The nonlinear diamagnetic signals above Tc are also 
observed in the single-layer Bi2Sr2_yLayCu06 (Bi 2201) 
family. In this series, the transition temperature Tc is 
tuned by the La content y. The optimal Tc occurs at 
y '^0.44. Samples with y > 0.44 are UD, while those 



with y < 0.44 are OV. 





FIG. 4; (color online) Curves of magnetization Md vs. H in 
Bi 2201 measured in intense fields. Panel A displays results 
on UD Sample Bi 2201 (La content y = 0.7 and Tc = 12 K). 
Note that most of the results shown are at T above Tc = 12 
K. They approach zero at the same nominal field scale~42 T. 
Panel B plots curves of Md vs. H measured in another UD 
Bi 2201 (La content y = 0.6 and Tc = 18 K). Panels C and D 
display Md curves measured in optimally-doped Bi 2201 (La 
content y = 0.44 and Tc = 30 K) and OV Bi 2201 (La content 
y = 0.2 and Tc = 20 K). In all these Bi 2201 samples, the 
curves of vs. H remain strikingly nonlinear high above 
Tc- In each panel, the curve closest to Tc is shown in bold. A 
preliminary version of Panel A was published in Ref. [ij • 

Figure |4] displays the diamagnetic Md — H curves of 
the single-layer cuprate Bi 2201 samples in the UD region 
(Panel A and B), in the optimally-doped region (C), and 
in the OV region (D). Above Tc, the Md — H curves in 
Bi 2201 are also similar to those in LSCO shown in Fig. 
[5] and Fig. [31 except that the magnitudes of Md and 
field scales are slightly smaller in the former. Above 
Md attains a broad minimum at fields below 20 T, and 
then approaches zero at > 40 T. Like the curves for 
LSCO in Fig. [2] and Fig. [31 the curvature of the low- 
field Mid — H curves changes from negative to positive 
as T increases across Tc- As H is increased beyond 20 



FIG. 5: (color online) Comparison of high-field magnetization 
curves Md vs. H in OPT Bi 2212 and UD Bi 2212. Panel A 
[13] displays the Ma curves in OPT Bi 2212 {Tc = 88 K) in 
fields up to 45 T. Panel B plots magnetization curves in UD 
Bi 2212 {Tc = 50 K) in fields up to 33 T. In each panel, the 
curve closest to Tc is shown in bold. Panel (A) is taken from 
ReL [H. 

T, Md is greatly suppressed. Above Tc, the complete 
suppression of Md requires very high fields - comparable 
to those needed below Tc. Even in UD Bi 2201 with y = 
0.7, where Tc = is quite low (12 K), the Md curves are 
suppressed to zero at ~ 38 - 42 T at T < 45 K. Above 
T^nset {"^ 50 K in this sample), Md vanishes throughout 
our entire field range. The interesting weak-field region 
is discussed under "fragile London rigidity" . 

In Fig. [51 we compare the curves in OPT and UD 
bilayer Bi 2212 (Panels A and B, respectively). Relative 
to the single-layer Bi 2201, the amplitude \Md\ in Bi 
2212 attains much larger values, and extend to higher 
field scales. By extrapolating the low-T Md-H curves, 
we estimate that Hc2 exceeds ^150 T (compared with 
50-80 T for Bi 2201). 

Fragile London Rigidity 

One of the most interesting features of the vortex liquid 
state above Tc is the fragile London rigidity, observable 
in the limit 7J — » 0. In Ref. [ll[, Li et al. discovered that 
over a broad interval of T (86-105 K) in OPT Bi 2212, 
the \ow-H Md follows the power-law dependence 

Md{T, H) ijiZ-siT) (ij ^ 0), (2) 

with an exponent 5{T) that grows rapidly from 1 (at 
T ~105 K) to large values (>6) as T ^ T+. This 
implies that the weak-field diamagnetic susceptibility 
X ~ limjj^o M/H — > — oo is weakly divergent through- 
out the interval in T where 6 > 1. However, this diver- 
gence is extremely sensitive to field suppression. The 
fragile London rigidity seems to refiect the increasing 
tendency of the phase-disordered condensate to establish 
long-range superfluid response as T — > Tj*" . It has no ana- 



5 



log in bulk samples of low- Tc superconductors, but may 
exist in a finite T interval above the Kosterlitz-Thouless 
(KT) transition in 2D systems such as Moi-ajGe^; and 

Using a soft cantilever, we have observed a similar pat- 
tern of magnetization in OPT Bi 2201. As shown in Fig. 

the Md curves display increasingly strong curvature 
as H approaches zero from cither direction. As T de- 
creases from 38 K to Tc (30 K), the zcro-i? slope rises 
sharply to a vertical line (see expanded scale in Panel 
B). The curve at Tc (bold curve) seems to approach a 
logarithmic dependence vs. H (equivalent to (5 ^ oo). 
(As may be seen by the oscillations, mechanical noise 
in this soft cantilever precludes accurate measurements 
for \H\ < 300 Oe. In Ref. [ul, high-resolution SQUID 
magnetometry was used to extend measurements down 
to 10 Oe, but the volume of the present Bi 2201 crystal 
is too small for similar SQUID measurements.) Despite 
the lower resolution, the divergent curvature apparent in 
Fig. [6] is consistent with the appearance of fragile Lon- 
don rigidity starting 8 K above T,.. The curves in Fig . 
[SJ3 are remarkably similar to those reported in [ll|, [3| 
for OPT Bi 2212. 

The fragile London rigidity is likely to extend over a 
larger T interval in UD samples. However, it would be 
more difficult to disentangle the intrinsic weak-if, non- 
linear M^-H behavior from the effects of inhomogeneous 
broadening arising from local variations of Tc. In OPT 
samples, we reason that such effects are minimized. 
The bulk of the sample has the maximum (OPT) Tc- 
Minority regions with lower Tc contribute only negligibly 
to the screening current. Hence, observation of the 
fragile state in OPT samples uncovers, in our opinion, a 
highly unusual feature of cuprates that is intrinsic. 
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FIG. 6: (color online) The low-field magnetization curves of 
OPT Bi2201 (La content y = 0.44). In Panel (A), the striking 
weak-field non-linearity is highlighted by displaying the varia- 
tion of Md from H= -2 T to +2 T. Although = 30 K (bold 
curve), diamagnetism is observed up to the onset tempera- 
ture T^nset —70 K. Panel (B) shows the low-field curvature in 
expanded scale. The measurements were performed using a 
very soft cantilever beam. The divergent curvature of Md at 
zero H is consistent with fragile London rigidity [HI . [l^ . 



of the magnitude \Md{T, H)\ is much larger as expected. 
Md{T, H) also extends to much higher field scales. A 
particular feature is that, near Tc, the contours in OPT 
Bi 2212 are nearly vertical up to 33 T (the maximum 
applied field) [l^. This implies that, at the separatrix 
at Tg, Md remains nearly H- independent up to 33 T. By 
contrast, in Bi 2201, the constancy extends only to 10 T, 
as may be seen in Fig. [71 



Contour Plot 

An instructive way to view the nonlinear diamagnetic 
magnetization is the contour plot of Md in the T-H 
plane [13]. Figure [7] displays the contour plot in single- 
layer UD Bi2201 (La content y = 0.7, Tc = 12 K). The 
value of \Md \ is as indicated at selected contours. With H 
fixed (e.g. at 10 T), \Md\ decreases monotonically as T is 
raised from 4 K to 60 K. Just as in the Nernst signal, the 
diamagnetic signal in the T-H plane bulges out to tem- 
peratures high above T^ with no obvious discontinuities 
or changes of slope. The highest temperature at which 
Md is resolved is '^50 K (the onset temperature in this 
sample). The absence of a boundary at Tc implies that 
the vortex-liquid state below Tc evolves continuously to 
the diamagnetic state above Tc- 

If we fix T (e.g. at 4 K) and increase 77, \Md\ also de- 
creases rapidly, as plotted in Fig. [H The field Hc2 (=45 
T) at which Md — > is plotted as open circles. With in- 
creasing r, HciiT) gradually decreases, roughly tracking 
the contour at 2 A/m. However, unlike the mean-field 
BCS scenario, Hc2 remains very large at Tc (arrow). The 
magnetization contour plots in Bi 2212 (see Ref. (Tsj ) 
are roughly similar to that in Fig. [T] However, the scale 



V. OPTIMALLY-DOPED YBCO 

Optimally-doped YBCO {Tc= 92 K) is distinguished 
as the cuprate with the smallest resistivity anisotropy 
and the largest interlayer (c-axis) coupling energy. Be- 
cause the coherence-length anisotropy S,a/(,c= 3-5 is only 
moderate, the vortices have the largest stiffness modu- 
lus along c among cuprates (^a and are the coher- 
ence lengths along the axes a and c, respectively). Ac- 
cordingly, the vortex solid melting line i7,„ {T) rises very 
rapidly below Tc (to ~15 T at 87 K). In the vortex-solid 
state [H < Hm), the dissipationless state survives to 
fields of 60 T or more. Optimally-doped YBCO should 
be the least susceptible to the phase disordering mech- 
anism for the destruction of long-range phase coherence 
at Tc (and hence the best candidate for Gaussian fluctu- 
ations among cuprates). 

However, the torque measurements reveal that Tc in 
OPT YBCO is also dictated by large phase fluctuations. 
Figure ^ displays the Md-H curves in OPT YBCO 
(twinned) measured to 14 T. The curves are broadly sim- 
ilar to those in LSCO 17 and OPT Bi2212, except for the 
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FIG. 7: (color online) Contour plot of the diamagnetic mag- 
netization \Md{T,H)\ of UD Bi 2201, with La content y = 
0.7 and Tc = 12 K (arrow). The spacing between adjacent 
contour lines is 10 A/m for T < Tc- The upper critical field 
Hc2 (defined by extrapolating Aid 0) is plotted as open 
circles. 



larger magnitude of \Md\ (at comparable H and T). At 
Tc= 92 K, \Md\ reaches the substantial value '^60 A/m 
at 14 T (by contrast, it should be nearly unobservable in 
a Gaussian mean- field picture). Over the broad interval 
92— s-130 K, a large diamagnetic signal is easily observed. 
As in the other hole-doped cuprates, very intense fields 
are needed to suppress in this interval. 

In Panel (B), we plot the T dependence of the total 
torque magnetization M^ff observed at 14 T (solid cir- 
cles). As in the other hole-doped cuprates, M^ff is un- 
resolved from the Van Vleck line {A + BT)H until T 
reaches T^^^^ ^130 K (arrow), below which it acceler- 
ates to very large negative values. For comparison, we 
have also plotted the magnetization measured in a very 
weak ~ 10 Oe (dashed curve) using a SQUID mag- 
netometer. The nearly vertical decrease signals flux ex- 
pulsion at Tc- We remark that, although Md - measured 
with H=10 Oe - is virtually unresolvable above Tc, the 
diamagnetic susceptibility x is actually quite large above 
Tc (as is clear from Panel A). Because is robust to 
intense H (100 T), the curve at 14 T reveals the exis- 
tence of the large fluctuating diamagnetism associated 
with the vortex liquid. This point, emphasized in Refs. 
[13, [HI, highlights the major difference between the dia- 
magnetism in hole-doped cuprates and low- Tc supercon- 
ductors. In the latter, increasing H in the fluctuation 
regime above Tc rapidly squelches the (Gaussian) fluctu- 
ation signal altogether. The curves in Fig. [51A display- 
ing significant diamagnetism surviving to intense fields, 
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FIG. 8: (color online). Magnetization curves of OPT 
YBa2Cu307-i. Panel (A) displays curves of Md vs. H from 
T = 89 to 150 K. Some curves were taken (at slower sweep 
rate) only to 1 T. Above Tc = 92 K (bold curve), a sizeable 



diamagnetic signal persists to Tj" 



130 K. Panel (B) plots 



the T dependence of the observed Me// with _ff = 14 T. As in 
the other cuprates, Afe// may be fit to a straight line above 
the onset of the diamagnetic signal Tj^^set- The dashed curve 
is the Meissner signal Msq measured by SQUID magnetome- 
try with H= 10 Oe. 



at temperatures up to 40 K above Tc is strong evidence 
that we are observing the phase-disordering mechanism, 
rather than Gaussian mean-field fluctuations. A com- 
parison of the Nernst and magnetization signals in UD 
YBCO is given in Ref. [H. 
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FIG. 9: (color online) Plots of the temperature dependence 
of M^ff{T) in Bi 2201 (Panel A) and in LSCO (B), showing 
the onset of diamagnetism as T is decreased. In both panels, 
the value of A/e// measured at ii" = 14 T is plotted vs. T in 
samples with various doping levels x. In general, Me// at high 
T varies weakly vs. T, as shown by the straight lines which 
are of the form A + BT. Relative to this linear background, 
Me// shows a strong downwards deviation starting at the 
onset temperature T^nset (indicated by arrows). 



VI. ONSET TEMPERATURES AND PHASE 
DIAGRAM 

An important question is how high in temperature does 
the diamagnetic signal extend above Tc- Following the 
procedure in Wang et al. [l^, we have plotted the total 
effective magnetization Meff measured in fixed H (14 
T) versus T. Figure [9] displays these plots for several 
samples of Bi 2201 (Panel A) and LSCO (Panel B). In 
agreement with the results for Bi 2212 in Ref. [l3l, Meff 
displays a weak T dependence at elevated T that may be 
fitted to the Van Vleck anisotropy term Axp = A + BT 
(straight lines). The anisotropy Axp leads to a paramag- 
netic torque background, as displayed. In each sample, 
Mf,ff(T) displays a sharp downward deviation, beginning 
at the temperature defined as T^nset (arrows). This re- 
flects the appearance of local supercurrents induced in 
response to the applied field. While feeble near T^g^n 
the diamagnetic term Md{T) grows very rapidly in mag- 
nitude to pull Meff to large negative values. 

To determine T^^^^t with reasonable accuracy, it is im- 
portant to establish the Van Vleck term Xp with a suf- 
ficiently dense set of points above T^nset- elevated 
temperatures, the magnetization curves become closely 
spaced. Thus, the stability of the cantilever as well as 
the resolution in measuring M^ff must be sufficiently 
high to allow adjacent curves to be distinguished. As ex- 
amples, we display in Fig. [TO] A and B expanded views of 
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FIG. 10: (color online) Expanded view of the total observed 
magnetization Me// vs. H in LSCO 09 (Panel A) and LSCO 
12 (B). Above T^nset, the T dependence of the curves is only 
from the paramagnetic Van Vleck term. Below T^set, how- 
ever, the diamagnetic term Md grows rapidly to dominate the 
T dependence. The stability of the torque cantilever and the 
resolution in t are sufficient to allow closely spaced curves to 
be resolved. The uncertainty in measuring Me// makes the 
largest contribution to the error bars in T^Let- 



the curves in LSCO 09 and LSCO 12, respectively. The 
values of M^f / at = 14 T are the ones plotted in Fig. 

m 

The profile of Meff vs. T is common to all the sam- 
ples investigated, even those in the extreme UD regime. 
The characteristic profile is qualitatively different from 
that seen in the gaussian regime in BCS superconduc- 
tors. Remarkably, the rapid downward acceleration of 
the diamagnetic signal matches the equally rapid growth 
of the Nernst signal taken at 14 T (Fig. 3 in Ref. [l3| 
compares the profiles of the Nernst signal and Md mea- 
sured in the same crystal of Bi 2212). 

To compare T^^set obtained here with the onset tem- 
perature of the vortex Nernst signal T^nset [2] , we plot the 
2 onset temperatures vs. doping x in the phase diagram 
for LSCO (Fig. [TTK) and Bi 2201 (Fig. [TTp). Remark- 
ably, in LSCO, T^nset (open squares) is nominally equal 
to T^nset ill LSCO over the entire doping range investi- 
gated. The major difference is that the former seems to 
peak at x = 0.12 whereas the latter peaks at 0.10. The 
error bars at both temperatures are too large to deter- 
mine if the disagreement is real. In the interesting UD 
side, both temperatures decrease roughly linearly with x 
as a; — > 0. [In an earlier analysis P], T^nset was extrapo- 
lated below X — 0.05 to reach at a; = 0.03. Our recent 
results show that this extrapolation is incorrect. Because 
the Nernst signal at a; = 0.03 is too weak to resolve even 
at low T, there is actually no experimental information 
on Tl^nsef By contrast, the results on Meff vs. T (Fig. 
[5J3) allow T^nset to be fixed reliably at small a;.] Inter- 
estingly, the X dependence of Hc2 obtained in Ref. [lB| 
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FIG. 11: (color online) Phase diagram comparing the on- 
set temperatures for the Nernst and diamagnetism sig- 
nals vs. doping X in La2-xSrxCu04 (Panel A) and in 
Bi2Sr2-yLayCu06 (Panel B). The superconducting transi- 
tion temperature Tc (sohd circles) is plotted with the on- 
set temperature Tonset determined by the Nernst effect (solid 
diamonds), and T^^^^ determined by torque magnetometry 
(open squares). In Panel B for Bi2Sr2-yLayCu06, a large 
La content y implies small hole carrier concentration (UD 
regime) . 



is also linear in x in this regime. 

In Bi 2201 (Panel B), the 2 temperature scales are also 
quite similar. However, the trend of T^^^s^t on the OV 
side appears to be slightly steeper than that of Tonset ■ ^ 
caveat is that the torque measurements here were not 
performed on the same crystals as the Nernst experi- 
ments. 

For the equivalent phase diagram of Bi 2212, see Ref. 
[13 . The phase diagram for YBCO appears in Ref. [24j . 



VII. DISCUSSION 

Diamagnetism and Supercurrent Response 
The presence of a large diamagnetic response that is both 
strongly T dependent and non-linear in H deeply impli- 
cates Cooper pairing. Diamagnetism involves an orbital 
current density J that is antiparallel to the applied vec- 
tor potential A (as in the London equation). Cooper 
pairing is - to our knowledge - the only established elec- 
tronic state capable of generating the current response 
consistent with the nonlinear, strongly T-dependent dia- 
magnetism reported here. [Core diamagnetism in insula- 
tors and Landau diamagnetism (observed in pure Bi) are 
both strictly H-hnear to extremely large H {hbH ~ W, 
where W is the band- width and /is the Bohr magne- 
ton) and nearly T independent. A "superdiamagnetic" 
state based on toroidal, orbital moments has been theo- 
rized ^T|, but this state has never been observed.] Hence, 
diamagnetism provides a rather direct detector of incip- 



ient Cooper pairing in the cuprates. 

As shown in Sees. HH JV] and |V] for LSCO, the 
Bi-based cuprates and YBCO, respectively, the Md-H 
curves above Tc show similar patterns - a broad mini- 
mum in moderate fields followed by a steady suppression 
to zero in very high fields. This pattern evolves continu- 
ously from the curves measured below Tc, which display 
a divergence at low field caused by the Meissner effect. 
Above Tc, this divergence vanishes because of the loss of 
long range phase coherence. Nonetheless, a reduced local 
phase rigidity survives [Tl| . which gives rise to the en- 
hanced diamagnetic Md above at low fields. Although 
its overall magnitude is ~10 times smaller than below 
(when observed at similar H), Md is readily detected as a 
strongly T-dependent and iJ-nonlinear response. In UD 
Bi 2201, we have accomplished full field suppression of 
Md in fields ~45 T. However, in all other cuprates, the 
full suppression requires fields in excess of 80 T (possibly 
as high as 150 T). These impressively large field scales are 
a consequence of the anomalously large binding energies 
of Cooper pairs in hole-doped cuprates. The broad sim- 
ilarity of the magnetization curves in LSCO, Bi 2201, Bi 
2212 and YBCO suggests that the diamagnetic behavior 
above Tc is universal in the hole-doped cuprates. They 
are qualitatively different from the diamagnetic response 
in low- Tc supercondcutors. 

It is instructive to compare the diamagnetism in 
cuprates with the fluctuating diamagnetism observed in 
disordered MgB2. In a recent experiment, Bernardi et 
al. 22] compared the magnetization of pure MgB2 (Tc = 
39 K) with disordered Mgi.^BaA^ (x = 0.25, Tc = 25 
K). In the disordered sample (which has a broad tran- 
sition width of ~15 K), the curves of Md vs. H show 
that sizeable diamagnetism exists in the narrow interval 
28-32 K above its Tc (Fig. 5 of [H). The profile of Md 
vs. H , which displays a broad minimum at ^ 200 Oe, is 
roughly similar to the profiles reported here (aside from 
the field scale) . The broad transition width implies large 
inhomogneities in the Al distribution. As the diamag- 
netic response above Tc does not persist above Tc of pure 
MgB2 (39 K), we suggest that the fluctuation diamag- 
netism arises from Al-poor regions of MgB2 which have 
the highest Tc. Thus, over the whole sample, diamag- 
netism is observable above 25 K, but not above 39 K. 
This contrasts with the cuprates. In the OPT sample 
within each family, Tonset extends above Tc by factors 
of 1.3 (YBCO), 1.4 (Bi 2212), 2.1 (LSCO) and 2.5 (Bi 
2201). Clearly, we cannot simply explain away the high- 
T diamagnetism as coming from isolated OPT regions 
with the highest Tc. The comparison shows that the lo- 
cal supercurrents detected in MgB2 arise from isolated 
regions with strong amplitude fluctuations and a broad 
distribution of local Tc's. By contrast, the diamagnetic 
signal in cuprates arises from a condensate that has lost 
phase stiffness, even though the gap amplitude remains 
large and nominally uniform above Tc. 

Quasiparticle term in Nernst signal 

In UD LSCO, the quasiparticle (qp) current makes a sig- 
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nificant contribution to the Nernst signal ■ In the ini- 
tial report of Xu et al. the onset temperature 
was found to remain high even when x falls below 0.1 
(T^ -150 K for a; = 0.05). This was traced ^ to a sig- 
nificant qp contribution to the Nernst signal. To discuss 
the qp term, it is crucial to consider the sign of the Nernst 
effect. 

By convention, the sign of the Nernst effect is defined 
as that of the triple product Ejv -H x (— VT) with E^v the 
observed Nernst i?-field [5| . This rule is equivalent to the 
old convention based on "Amperian" current direction 
(clearly described by Bridgman For vortices, E^v = 

B X Vi where the vortex line velocity vj;, is ||(— VT). 
Hence vortex flow produces a positive Nernst signal. The 
qp contribution may have either sign (unrelated to their 
charge sign). 

When the qp term is negative, it is relatively easy to 
separate the 2 contributions, especially by going to in- 
tense H . However, if the qp term is positive, the separa- 
tion is more difficult. In Ref. Wang et al. introduced 
a method for separating the qp and vortex terms by si- 
multaneous measurements of the Hall angle Oh and ther- 
mopower S vs. H to obtain the term StanOH- This sub- 
traction procedure yields the onset temperature T^nset 
for the vortex term (45 K at x = 0.05), which is plot- 
ted in Fig. [TTK. With the qp subtraction applied to the 
samples x < 0.07, the curve of T^nset ^ has a tilted 
dome profile with a sharp peak at x = 0.125. For OPT 
and OV LSCO, Wang et al. found [5] that the qp term is 
negative and negligible compared with the vortex term. 
Hence, T^nset gives the onset of the vortex term without 
the need for corrections. (The qp term is also negative 
in YBCO, Bi 2201 and Bi 2212 [5]). 

Recently, this assumption has been challenged by Cyr- 
Choiniere et al. [18], who proposed that rearrangement of 
the FS - possibly by charge ordering or stripe formation 
- produces a positive quasi-particle term that dominates 
the Nernst signal over the LSCO phase diagram at high 
temperatures. Although their Nernst measurements were 
largely on Eu-doped and Nd-doped LSCO where static 
stripes are experimentally observed, they have extended 
their hypothesis to pure LSCO. There, they proposed 
that a positive qp also accounts for its Nernst signal at 
elevated T (except in a narrow interval just above Tc). 

In the scenario of Ref. [l^ , the onset temperature for 
vortex fluctuations should lie considerably lower than 
the dome of T^nset plotted in Fig. [TTK. However, the 
good agreement between T^^set ^^'^ "^onset shows that 
this is not the case. As we argued above, the strongly 
T- and if-dependent diamagnetism arises only from the 
pair condensate, and is unaffected by qp contributions. 
The agreement between T^^^^ and T^nset seems to us to 
be strong evidence against the proposal in Ref. [l^, at 
least in pure LSCO. For the claim to be viable, the hypo- 
thetical quasiparticles would have to produce a large, T- 
dependent diamagnetism that is also strongly non-linear 
in _ff , as well as a positive Nernst signal. 

The hypothesis of Cyr-Choiniere et al. [l^l is the lat- 



est of several proposals (see discussion in Ref. ^) that 
have sought to explain away the unexpected Nernst sig- 
nal in the cuprates by invoking quasiparticles with ad hoc 
properties. Are the large Nernst signals at high T from 
vortices in a phase -disordered condensate or from quasi- 
particles (qp)? We argue that torque magnetometry and 
the Nernst effect together constitute an incisive combi- 
nation that answers this question. When the 2 probes 
show that a large positive Nernst signal coexists with a 
diamagnetic susceptibility (with the same onset temper- 
atures and similar proflles vs. T and H), the case in favor 
of phase slippage in a pair condensate with strongly dis- 
ordered phase seems compelling to us. This is one of our 
main conclusions. 

The magnetization-Nernst approach can also be turned 
around to identify situations when the Nernst signal is 
not caused by vorticity. As mentioned, in UD LSCO 
(x <0.07) the agreement of T^^^^ and T^nset provides 
confirmation that the subtraction proceduce based on 
S tan 9h is valid. (I n the broad interval between T^^set 
and T,y, extending from 45 to 150 K for a; = 0.05, there is a 
large positive Nernst signal, but diamagnetism is absent. 
Thus, even if the qp contribution had not been identi- 
fied [2], the present experiment would have detected the 
correct onset of the vortex term.) 

In UD YBasCusOg+y (YBCO), the qp Nernst signal 
is unusually large in a narrow window of doping. In this 
doping range, when the qp term appears at high T, it is 
negative as in OV LSCO. Signiflcantly, the diamagnetic 
signal is absent until the vortex Nernst signal appears 
at a lower T (detailed YBCO results are reported in 
Ref. [13] )■ Hence, when used together, the torque 
magnetization and Nernst effect readily distinguish qp 
from vortex contributions to the Nernst signal. Nernst 
effect and diamagnetism studies on La2-xBaxCu04 (20| 
should provide valuable insight into the fluctuation 
regime in which stripes coexist with Cooper pairing. 

Related Experiments 

The surviving pair condensate above Tc has also been ob- 
served in other experiments, notably in measurements of 
the kinetic inductance [1^ , STM experiments on the gap 
above Tc |26| . and survival of Bogolyubov quasiparticles 
above Tc [23]. Consistent with these observations, our 
results imply that the pair condensate exists well above 
Tc, surviving as a dilute vortex liquid with local phase 
rigidity of short phase-correlation length. A number of 
groups recently calculated the Nernst signal and diamag- 
netism above Tc in 2D superconductors and applied the 
results to cuprates [.28, 29, 30, 31, 32]. The vortex hquid 
viewed as an incompressible superfluid has been treated 
by Anderson [H, [sj] • A relevant discussion of the rela- 
tion of the quantum oscillation results to the Nernst and 
magnetization results is given in Ref. [s^] . 

To summarize, the high-fleld torque magnetometry 
measurements reveal that the diamagnetism persists well 
above Tc in several families of hole-doped cuprates. As 
a strongly nonlinear Md vs. H is characteristic of local 
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supercurrent response, the diamagnetism is direct evi- 
dence that the pair condensate exists above Tj,, surviv- 
ing in places all the way to T^nset- Hence, phase slippage 
is the origin of the large Nernst effect signals observed 
to that temperature. The agreement between the onset 
temperatures T^nset ^^'^ "^onset precludes quasiparticle 
interpretations for the positive Nernst signal above Tc- 
The magnetization results pose very serious difficulties 
for the quasiparticle hypothesis [l^l proposed recently 
for the Nernst effect in pure LSCO. 

We acknowledge numerous helpful discussions with P. 



W. Anderson, J. C. Davis, S. A. Kivelson, P. A. Lee, 
T. Senthil, Z. Tesanovic and A. Yazdani. The research 
at Princeton is supported by funds from U.S. National 
Science Foundation under the MRSEC Grants DMR- 
0213706 and DMR-0819860. Y.W. is supported by NSFC 
and MOST of China. G. D. G. is supported by the 
Department of Energy (DOE) under contract No. DE- 
AC02-98CH10886. The high-field experiments were per- 
formed at the National High Magnetic Field Laboratory, 
which is supported by NSF Cooperative Agreement No. 
DMR-084173, by the State of Florida, and by the DOE. 



Z. Xu, N.P. Ong, Y. Wang, T. Kakeshita and S. Ucliida, 
Nature 406, 486 (2000). [17] 
Yayu Wang, Z. A. Xu, T. Kakeshita, S. Ucliida, S. Ono, 
Yoichi Ando, and N. P. Ong, Pliys. Rev. B 64, 224519 [18 
(2001). [19 
Yayu Wang, N. P. Ong, Z.A. Xu, T. Kakesliita, S. 
Uchida, D. A. Bonn, R. Liang and W. N. Hardy, Phys. 
Rev. Lett. 88, 257003 (2002). ' [20 

Yayu Wang, S. Ono, Y. Onose, G. Gu, Yoichi Ando, 
Y. Tokura, S. Uchida, and N. P. Ong, Science 299, 86 [21 
(2003). 

Yayu Wang, Lu Li, and N. P. Ong, Pliys. Rev. B 73, [22 
024510 (2006). 

V. J. Emery, S. A. Kivelson, Nature 374, 434 (1995). [23 

C. Capan, K. Belmia, J. Hinderer, A. G. M. Jansen, W. [24 
Lang, C. Marcenat, C. Marin, and J. Flouquet, Phys. 

Rev. Lett. 88, 056601 (2002). [25 
F. Rulfier-Afbenque, R. Tourbot, H. Alloul, P. Lejay, 

D. Cofson, and A. Forget, Phys. Rev. Lett. 96, 067002 [26 
(2006) 

F. Rulfier-Afbenque, H. Affouf, Cyrif Proust, P. Lejay, 
A. Forget, and D. Cofson, Phys. Rev. Lett. 99, 027003 [27 
(2007). 

Yayu Wang, Lu Li, M. J. Naughton, G. Gu, S. Uchida 
and N. P. Ong, Phys. Rev. Lett. 95, 247002 (2005). [28 
Lu Li, Yayu Wang, M. J. Naughton, S. Ono, Yoichi Ando, 
and N. P Ong, Europhys. Lett. 72, 451 (2005). 
Lu Li, Yayu Wang, J.G. Checkefsky, M.J. Naughton, [29 
Seiki Komiya, Shimpei Ono, Yoichi Ando, and N.P. Ong, 
Physica C 460-462, 49 (2007). [30 
Lu Li, Yayu Wang, M.J. Naughton, Seiki Komiya, Shim- [31 
pei Ono, Yoichi Ando, N.P. Ong, J. Magnetism Magnetic 
Mater. 310, 460 (2007). [32 
N. P. Ong, Yayu Wang, Lu Li, and M. J. Naughton, Phys. 
Rev. Lett. 98, 119702 (2007). [33 
Lu Li, J. G. Checkefsky, Seiki Komiya, Yoichi Ando and [34 
N. P. Ong, Nature Physics 3, 311 (2007). [35 
D. C. Johnston, and J. H. Cho, Phys. Rev. B 42, 8710 



(1990). 

A. Lasciaffari, A. Rigamonti, L. Romano', A. A. Var- 
famov, and I. Zucca, Phys. Rev. B 68, 100505(R) (2003). 
Ofivier Cyr-Choiniere et al, Nature, 458, 743 (2009). 
C. Bergemann, A. W. Tyfer, A. P. Mackenzie, J. R. 
Cooper, S. R. Jufian, and D. E. Farreff, Phys. Rev. B 
57, 14387 (1998). 

Q. Li, M. Hucker, G. D. Gu, A. M. Tsvelik, and J. M. 

Tranquada, Phys. Rev. Lett. 99, 067001 (2007). 

V. L. Ginzburg, A. A. Gorbatsevich, Yu. V. Kopayev and 

B. A. Vofkov, Sofid State Commun. 50, 339 (1984). 

E. Bernardi, A. Lasciaffari, A. Rigamonti, and L. Ro- 
mano', Phys. Rev. B 77, 064502 (2008). 
P. W. Bridgman, Phys. Rev. 24, 644 (1924). 
Minhao Liu, Yayu Wang, Lu Li and N. P. Ong, in prepa- 
ration. 

J. Corson, R. Maffozzi, J. Orenstein, J. N. Eckstein, and 
L Bozovic, Nature, 398, 221 (1999). 
Kenjiro K. Gomes, Abhay N. Pasupathy, Aakash Pushp, 
Shimpei Ono, Yoichi Ando, Afi Yazdani, Nature, 447, 
569 (2007). 

Jhinhwan Lee, K. Fujita, A. R. Schmidt, Chung Koo 
Kim, H. Eisaki, S. Uchida, and J.G. Davis, Science 325, 
1099 (2009). 

L. Benfatto, C. Casteffani, T. Giamarchi, Phys. Rev. 
Lett. 98 117008 (2007); ibid., Phys. Rev. Lett. 99, 207002 
(2009). 

V. Oganesyan, D. A. Huse, and S. L. Sondhi, Phys. Rev. 
B, 73, 094508 (2006). 

Z. Tesanovic, Nature Physics 5, 408 (2008). 

Danief Podofsky, Srinivas Raghu, and Ashvin Vish- 

wanath, Phys. Rev. Lett. 99, 117004 (2007). 

Subroto Mukerjee and David A. Huse, Phys. Rev. B 70, 

014506 (2004). 

P. W. Anderson, Nature Physics 3, 160 (2007). 

P. W. Anderson, Phys. Rev. Lett. 100, 215301 (2008). 

T. Senthif, P. A. Lee, Phys. Rev. B 79, 245116 (2009). 



